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Mechanical exfoliation is a convenient and effective approach to deriving two-dimensional (2D) 
nanodevices from layered materials; but it is also generally perceived as unpreferred as it often yields devices 
with structural irregularities and nonidealities. Here we show that such nonidealities can lead to new and 
engineerable features that should be embraced and exploited. We measure and analyze high frequency 
nanomechanical resonators based on exfoliated 2D molybdenum disulfide (MoS 2 ) structures, and focus on 
investigating the effects of structural nonidealities and asymmetries on device characteristics and 
performance. In high and very high frequency (HF/VHF) vibrating MoS 2 devices based on diaphragms of 
—2—5 fim in size, structural nonidealities in shape, boundary, and geometric symmetry all appear not to 
compromise device performance, but lead to robust devices exhibiting new multimode resonances with 
characteristics that are inaccessible in their 'ideal' counterparts. These results reveal that the seemingly 
irregular and nonideal 2D structures can be exploited and engineered for new designs and functions. 



Two-dimensional (2D) nanostructures in graphene and transition metal dichalcogenides (TMDCs), derived 
from their corresponding layered materials in bulk, have recently attracted significant attention due to their 
unique properties for creating new devices with potential for nanoelectronics, optoelectronics, and trans- 
ducers 1 " 4 . In device fabrication, it was mechanical exfoliation that helped make a breakthrough in isolating 
atomically thin 2D devices from layered materials, yielding graphene from graphite 5 . The seemingly simple 
mechanical exfoliation method comes along with great intuition and insight - it takes advantage of the fact that 
atomic sheets/planes are weakly bonded by van der Waals interactions in layered materials such as graphite and 
TMDCs. In practice, its convenience, cleanness (especially in the device region), and low cost (e.g., via Scotch 
tapes) add to its merits. Following graphene, atomically thin hexagonal boron nitride (h-BN), molybdenum 
disulfide (MoS 2 ) and other TMDCs have been successfully obtained using exfoliation, in realizing the first 2D 
devices from these crystals 6 . Today, 'handcraft' of mechanical exfoliation is still an important and effective 
approach to yielding pristine and high-quality 2D nanodevices, including sophisticated heterostructures 7 " 11 . 
While exfoliation can be limited by hand operations in practice (perhaps primarily on being at small scale) today, 
there is great room in advancing scalable exfoliation-type machinery processes to systematically and fully exploit 
the weak van der Waals interactions in crystals of layered materials. 

In current prevailing 2D nanostructures and devices made by exfoliation methods, structural nonidealities in 
shape, boundary, and geometric asymmetry are known to exist. These nonidealities have hitherto usually been 
deemed undesirable; and such nonideal-looking structures are generally discarded. There has been no systematic 
study on how the nonideal appearances and geometric irregularities actually affect these structures as devices. 

In this work, we use nanomechanical resonators made of ultrathin suspended MoS 2 layers as prototypical 2D 
devices, and systematically investigate the effects of structural nonidealities on their attributes and performance. 
We choose MoS 2 because it is a 2D semiconductor (different than graphene) with band structure depending on 
number of layers (thickness) and tunable via strain 12 " 14 . Lately it is also suggested that MoS 2 resonators may be 
intrinsically less dissipative than graphene counterparts owing to reduced phonon-phonon damping 15 . The 
nonidealities we focus on include irregular shapes/edges and asymmetric geometries. From extensive measure- 
ments and modeling, we demonstrate that "nonideal" MoS 2 structures can actually make robust and uniquely 
interesting multimode nanomechanical resonators, offering new opportunities for both fundamental studies and 
novel technological applications; we also explain the observed effects and develop guidelines for engineering such 
devices toward potential utilization. 
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Suspended MoS 2 flakes are made by mechanical exfoliation of 
MoS 2 onto prefabricated microtrenches, to produce MoS 2 nanome- 
chanical resonators with a diaphragm- on-microtrench drumhead 
structure (see Methods and Supplementary Information SI). 
Devices with different types of geometrical nonidealities can be 
obtained. To directly measure these pristine exfoliated structures 
and avoid possible contaminations in any extra processing steps 
{e.g., coating with polymer, patterning contact electrodes), here we 
employ laser interferometric techniques that were often used for 
nanomechanical resonators in conventional materials 16 " 20 ; and we 
have further carefully engineered the interferometric readout tech- 
niques (with significant improvement in displacement sensitivity) so 
that we are able to directly observe the undriven (no external drive), 
intrinsic thermomechanical (Brownian) motions of our nanoscale 
2D structures 21,22 (see Methods and Supplementary Information S2 
for details). 



Results and discussion 

Incomplete drumhead devices with free edges. We first investigate 
the effects due to the nonideality of incomplete MoS 2 coverage over 
the microtrenches. When a microtrench is not completely covered, 
we have a partial drumhead structure (Fig. la and b) whose perimeter 
is partly undamped, exhibiting an open arc and a free edge (a chord 
corresponding to the open arc). As illustrated in the SEM images, 0 is 
the central angle of the open arc. We call 0 the "open angle", which 
measures the degree of incomplete coverage of the microtrench. 
Despite its geometry of being an incomplete drumhead, such 
structures still make robust nanomechanical resonators. The 
resonances are found to fit well to the damped harmonic resonator 
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Supplementary Information S3 for details), with quality (Q) factors 
comparable to Q values of the devices based on their ideal, 
completely- covered counterparts (see Supplementary Information 
S4 & Supplementary Table SI for details of multiple devices). 

We analyze the effects of incomplete drumhead structure by using 
finite element modeling (FEM, in COMSOL). Figure lc shows the 
computed resonance frequencies and mode shapes for the 3 lowest 
flexural modes (fi — / 2a — »^>/2b — »"^)> as functions of the 
open angle 0, for a partially clamped device (with diameter d = 6 um, 
and thickness t = 30.3 nm). Note here/x is the fundamental mode 
with no nodal line in the mode shape, / 2a is the mode with a single 
nodal line in parallel to the free edge, and/ 2b is the mode with a single 
nodal line perpendicular to the free edge, respectively. A few inter- 
esting features are observed. First, the resonance frequency is not a 
monotonic function of open angle: as 0 increases,^ 0) first decreases, 
then increases. All the 3 lowest modes exhibit this pattern. Second, 
the shapes of thefi(O) curves, such as the open angle 0 m at which the 
frequency assumes its minimum f i>m = fj(6 m ), depends on the indi- 
vidual mode (with index t). This leads to mode crossing (f 2a values 
surpassing/ 2b value) at 0 ~ 137°. Third, the shape of each curve does 
not show strong dependence on the device dimensions. Shown in 
Fig. Id for a larger thin device (d = 10 urn, t = 10 nm) and a smaller 
thick device (d = 1 um, t = 100 nm), the solid and dashed curves 
overlap each other (after rescaling), despite the 100 times difference 
in aspect ratio d/t (diameter/thickness). These observations on the 
fi(6) curves open a new possibility of mode shape engineering by 
designing and tuning the open angle 0 of an incompletely clamped 
diaphragm. As demonstrated in Fig. lc, for fully covered devices (6 = 
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Figure 1 | Characteristics and features of nanomechanical MoS 2 resonators with incomplete drumhead geometries, (a): SEM image (top), optical 
image (bottom), and measured mechanical resonance of a device with diameter d = 2.53 (am, thickness t = 43.1 ± 1.4 nm, and open angle 6 ~ 100°. (b): 
Images and measured data from another device with d = 2.69 um, t = 62.2 ± 0.7 nm, and 6 ~ 102°. Scale bars (for SEM images only): 2 (am. (c): The 
computed f versus 6 curves (6 resolution: 2°) and mode shapes of the three lowest modes for a device with d = 6 (am and t = 30.3 nm (dimensions based 
on device #4 in Supplementary Information Table SI). The computed mode shapes are shown for selected 6 values (integer multiples of 30°, indicated by 
the open squares on the curves), (d): Comparison of the f— 6 curves between a small/thick device (solid curves, d = 1 (am, t = 100 nm, left y axis) and a 
large/thin device (dashed curves, d = 10 (am, t = 10 nm, right y axis), (e): Dependence of the clamping length ( oc fc eff ), resonator area ( oc M eff ), and 
resulting frequency / oc (/c eff /M ef f) 1/2 on the open angle 6 (all values normalized to the 6 = 0 case, i.e., the fully clamped case). 
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0 only), one is limited in the frequency ratios between the modes, and 
in degenerate mode pairs with small and random frequency splitting. 
With 0 7^ 0, one can adjust the 6 value and choose many possible 
frequency ratios between different modes, and the mode sequence 
can also be changed by engineering the 9 values. 

The observed fi(6) behavior can be understood, very intuitively, 
by examining the interplay between the changes in both clamping 
conditions and resonator size, as 9 increases. For a circular 
diaphragm clamped at the edge with flexural rigidity D, and 
deflection w{r,(p), the associated potential energy is 23 U = 

\dr 2 



im: 



IdwX 
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rdr for mode shapes that are azimuth- 



ally symmetric. By taking spatial derivatives we obtain an effective 



stiffness of k 



'eff-- 



d 2 

dw eff 2 



U-- 



dcp x 



D d 2 

2 dw eff 2 



d 2 w 



1 dw 



dr 2 r dr 



rdr, 



where w e ff is the effective displacement of the entire diaphragm. For a 
partially clamped diaphragm, by analogy we approximate the effec- 
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a central angle (2% — 0)) provides all the restoring force. At small 0 
values, increasing 0 mainly reduces the clamping arc length. As 0 
continues to increase, the reduction in device area (and thus the 
effective mass M e ff) becomes more dominant, and the frequency 
increases again (f oc [A; e ff/M e f f ] 1/2 ). Behavior of the higher modes 
can be understood similarly, with different /c eff and M eff values for 
different modes of the same device structure (see Supplementary 
Information S5 for detailed discussion). Figure le plots the value 
of [A: e ff/M e f f ] 1/2 as a function of 0. Clearly, between 0 and n, the 
frequency first decreases and then increases, capturing the main 
feature of this phenomenon. 

From Fig. Id one can see that the effect of open angle 0 upon the 
frequency curves is purely geometric - the shapes of fj{6) curves only 
depend on 0 and do not depend on actual dimensions of the devices 
(see Supplementary Information S6 for additional discussion). This 
fact allows one to quantitatively determine the device geometry (0 
value) from the frequency ratios alone, regardless of device size. 
Figure 2a shows a device with 0 ~ 160° (measured from SEM and 
optical images), whose measured resonance spectrum is shown in 
Fig. 2b. The calculated frequency ratios for 0° < 0 < 180° is plotted in 
Fig. 2c. By comparing calculated and measured frequency ratios we 



determine a 0 value in the range of 155° — 165°, in good agreement 
with the actual geometry of the device. 

The fact that such partial drumhead devices with free edges 
can have comparably high resonant performance as complete 
drumheads structures has important implications. First, this imme- 
diately increases the yield of functional devices by a factor of 2 or 
more (today in making suspended graphene and MoS 2 2D structures 
by mechanical exfoliation, often —50% suspended structures may 
have partially- covered geometries with free edges, and yielding sus- 
pended MoS 2 devices is particularly more challenging). More inter- 
estingly, edges of MoS 2 crystal have distinct chemical properties 
different than the 2D surface, and have been studied for applications 
such as electrocatalysis 24 28 . This enables one to preferentially func- 
tionalize the MoS 2 resonator and create binding sites at selective 
locations (e.g., along the edges), offering opportunities for loca- 
tion-specific chemical and biological sensing. Furthermore, the edges 
of 2D crystals have unique electrical properties which make them 
interesting ID sub-systems. For example, graphene edges have chir- 
ality- dependent electronic states 29,30 , and can remain gap-less even 
when a symmetry-breaking-caused energy gap is opened in the 2D 
bulk 31 . Such edge states are expected to carry dissipationless spin 
currents, as predicted from the topological properties of graphene 
as a quantum spin Hall insulator under spin-orbit coupling 31 " 34 . 
The properties of such spin-filtered ID channels are expected to 
be sensitive to perturbations such as flexural curvature 33 , which 
can naturally couple to resonant motions in these partially- covered 
drumhead structures. In MoS 2 , the metallic ID edge 35 can be 
ferromagnetic for certain atomic configurations 36 . Both graphene 
and MoS 2 edge states have been experimentally observed in 
substrate- supported samples 35,37 ' 38 . With ID electrical contacts 
directly accessing the edge of such 2D crystals 39 , these partial- 
drumhead resonators may soon help make new platforms for study- 
ing the coupling between mechanical, electrical, and spin degrees 
of freedom in these ID edge states, such as effects due to curvature 33 
and the interaction between spin-orbit and electron-phonon 
couplings 40 . 

Asymmetry enables new devices - elliptic versus circular. We now 

explore the effects of geometric asymmetry and shape distortion of 
devices, by examining elliptic versus circular devices. Such imperfec- 
tions may occur, sometime systematically, with challenging designs 
or imperfections in the manufacturing processes. Figure 3a shows 
three MoS 2 diaphragms with the same thickness but different shapes 
- a circle and two ellipses. All three devices exhibit clear well-defined 
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Figure 2 | Quantifying device geometry through measured frequency ratios of multimode resonances from incomplete drumhead devices, (a): SEM 
(perspective and top views) and optical images of a MoS 2 resonator with incomplete drumhead structure. From the images we measure 6 ~ 160°. (b): The 
measured thermomechanical noise spectrum of the device, showing the first three resonance modes. Insets: Zoom-in on the individual resonances 
(measured at respective optimized locations on the device) showing the fitting to the model, (c): The comparison between calculated (solid curves) and 
measured frequency ratios between different resonance modes (measured values indicated by the horizontal dashed lines). The highlighted area denotes 
the 6 range (155°-165°) in which the measured values intercept the calculated ones, showing good agreement with the 6 value determined from the 
images. 



SCIENTIFIC REPORTS | 4:3919 | DOI: 1 0.1 038/srep0391 9 



3 



b 1.0 

r0.8 




(D 0=250 


© £=210 




I u 


: 0 

















14.0 14.5 17.0 17.5 11.5 12.0 
Frequency (MHz) 




Frequency (MHz) 

Figure 3 | Elliptic MoS 2 nanomechanical resonators and mapping of their multimode resonances, (a): Optical and SEM images of three elliptic and 
circular resonators based on the same MoS 2 nanosheet (t = 57.5 nm) (b): Measured thermomechanical resonances from device 1 (a = 6.81 um, b = 
5.89 um), device 2 {d = 5.53 um), and device 3 (a = 6.78 (am, b = 5.89 (am), (c): Wide range frequency spectra measured from both ecliptic devices (1 
and 3). Insets: Spatially mapped resonance mode shapes for the first three resonance modes, measured by the scanning interferometric spectromicroscopy 
technique, in comparison to the simulated mode shapes in finite element modeling. Dashed lines outline the boundaries of the elliptic drumhead 
resonators. 



fundamental-mode resonances (Fig. 3b). The measured Qs show no 
observable dependence on the in-plane aspect ratio of the device (al 
b) while the resonance frequency clearly increases with decreasing 
device size. 

When measured over a wider frequency range, both elliptic 
devices exhibit multimode resonances (Fig. 3c). We further spatially 
map the mode shapes of the observed multimode resonances, by 
performing the 'scanning interferometric spectromicroscopy' mea- 
surements, in which we focus on each individual resonance and scan 
the laser spot over the entire device surface (with a sub -micron 
spatial resolution), and measure the resonance amplitude versus 
the laser spot position. Figure 3c insets show the measured mode 
shapes vis-a-vis FEM simulations for all the resonances, demonstrat- 
ing good agreement. We confirm that, due to the elliptic shape, the 
degeneracy in the 2 nd and 3 rd modes is lifted, and the mode with nodal 
line along the semi-major axis has higher frequency than the mode 
with nodal line along the semi-minor axis. 

To fully understand how the deviation from a circle to an ellipse 
affects the resonance frequencies, we perform FEM simulations for 
elliptic MoS 2 diaphragm resonators with different aspect ratios {alb, 
the ratio in length between the two semi- axes, see illustrations in 
Fig. 3b and Fig. 4). While resonance frequencies of elliptic diaphragms 
can be solved analytically with some approximation 41,42 , here we focus 
on the evolution of mode shapes and mode crossing. For ellipses, we 
use the notation pXq, the arrangement of anti-nodal motional 
domains (individual singly-connected continuous areas enclosed by 
nodal lines), to describe the mode shapes (see Fig. 4 legends). Note that 



this notation for ellipses is consistent with the conventions for describ- 
ing multimode vibrations in rectangular diaphragms 43 . For circular 
devices, we use the traditional notation (n,m), where n and m refers 
to the number of nodal diameters and nodal circles (Fig. 4 legends). 

From Fig. 4 we find that the lowest resonance frequency mode (dark 
yellow curve) is always the elliptic 1X1 mode, which corresponds to 
the circular (0,1) mode. The next lowest ones are the elliptic 1X2 and 
2X1 modes (red/blue curves), corresponding to the circular (1,1) 
mode. At alb = 1 these two modes degenerate and the curves cross 
each other. All these 3 modes are observed in Fig. 3c for both elliptic 
devices, showing good agreement between measurement and theory in 
both mode sequence and mode shapes. 

The next lowest modes near alb = 1, 1X3 and 3X1 (wine and 
magenta curves) show more complicated behavior. As alb 
approaches 1, their mode shapes deform, with the center motional 
domains elongating into dumbbells and eventually bifurcating, 
becoming essentially the elliptic 2X2 modes (olive curve) at alb = 
1 (corresponding to the circular (2,1) mode). In contrast to the 1X2 
and 2X1 modes, where the mode shapes continue on both sides of al 
b = 1, here the 1X3 and 3X1 modes on one side do not continue into 
the same mode shapes on the other side (double-headed arrows in 
Fig. 4). Indeed, the higher frequency 3X1 and 1X3 modes (magenta- 
wine curve), whose three motional domains align along the semi- 
minor axis, merge into the circular (0,2) mode at alb = 1. 

The evolution of the mode shapes also affects how the curves 
approach each other upon degeneracy. When the mode shapes are 
retained in both regimes across alb = 1, as in the 1X2-2X1 cross- 
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Figure 4 | Frequency scaling of the multimode resonances in the elliptic MoS 2 nanomechanical resonators. The frequency versus aspect ratio plot shows 
computed results for an elliptic MoS 2 resonator with t = 57.5 nm, b = 5.53 um (dimensions based on device (2) in Fig. 3), and alb varying from 1/2 to 2 
(with device shape illustrated for selected al lvalues under the horizontal axis). The device orientation is shown on the alb = 1.5 case under the axis and on 
the 1X4 mode shape at the top of the plot. Each curve (with distinct color and line type) shows the calculated /values for a particular mode shape 
(labeled on the corresponding curve). Double-headed arrows indicate the anti-crossing between the 1X3 and 3X 1 modes at alb = 1. Legend to the right 
gives the naming convention for the elliptic ( top) and circular ( bottom) mode shapes, with dashed lines representing the nodal lines. Inset At higher 
frequencies, the 1X4/4X1 modes (blue curve, with extrapolations shown in dotted lines) and the 2X3/3X2 modes (red curve) cross at alb = 1 with 
degeneracy in both frequency and mode shape. 



ing, the curves (red-blue) simply cross, retaining individual slopes. 
When the mode shapes also degenerate (one mode shape evolves into 
another) , as in the 1X3-2X2 and 1X4-2X3 (Fig. 4 inset) crossing, 
the curves become tangent to each other. When a certain mode 
evolves into different shapes when approaching alb = 1 from oppos- 
ite directions, such aslX3^2X2 from alb < 1 and 1 X 3^ (0,2) 
from alb > 1, anticrossing is observed (doubled-headed arrows in 
Fig. 4). 

By examining the mode evolutions, we clearly demonstrate and 
verify that circular devices are in fact one special case of elliptic 
devices, with some circular modes being the degenerate case of mul- 
tiple elliptic modes. As a result, any broken symmetry in shape on a 
circular resonator can lift its degeneracy and dictate the resonant 
modes to take particular orientations, split in frequency, or even 
bifurcate into different elliptic mode shapes. Such effects are prac- 
tically unavoidable in intentionally designed circular resonators, 
leading to random, undesirable splitting of the degenerate modes. 
The results here suggest that designing the resonators to be elliptic 
can alleviate such problems by lifting the degeneracy in a controlled 
manner. 

Lifting resonance degeneracy is important for realizing next-gen- 
eration devices for spatially- resolved resonant sensing of surface 
adsorbates. Study of surface adsorption on ID NEMS devices has 
already revealed a number of interesting phenomena such as low- 
dimensional phase transition 44 ' 45 and surface adsorbate fluctua- 
tions 46 ' 47 . By monitoring multiple resonance modes of a MEMS/ 
NEMS resonator, it is possible to simultaneously determine the mass 
and position of the adsorbed species 48 " 51 . Such functionality, so far 
only demonstrated in ID resonators, could be extended to 2D res- 
onant sensors 52 . Spatially- resolved detection and dynamical moni- 
toring of surface adsorbates appear to be more interesting on 2D 



platforms because the additional degree of freedom allows for richer 
position responsivity coupled to multiple 2D mode shapes, thus 
enabling studies of intriguing phenomena such as on -device cell 
division and 2D chemical reactions 53 with simultaneous mass and 
position detection of individual particles. In order to unambiguously 
determine the adsorbate position, it is desirable that any random 
degeneracy be removed from the multimode characteristics of the 
resonant sensor. In this case, perfectly circular resonators become 
unfriendly for resolving positions of adsorbates, for having the high- 
est degree of symmetry and degeneracy among 2D geometries. 
Elliptic resonators demonstrated in Fig. 4 systematically lift the 
undesired resonance degeneracy, and are thus better suited for multi- 
mode resonant sensing. We envisage, exploring device geometries 
with even higher degree of asymmetry (e.g., specially- engineered 
scalene triangular resonators) can further control the mode shapes 
of degenerate resonances, completely removing the spatial ambiguity 
in 2D multimode monitoring of adsorbates. 

Conclusions 

In summary, we have experimentally demonstrated that MoS 2 2D 
nanomechanical resonators with structural nonidealities and asym- 
metries, compared with their 'ideal' counterparts, can engender new 
multimode resonance characteristics. Given the occurrences in large 
populations of exfoliated devices made of emerging 2D materials 
(including graphene, h-BN, and various TMDCs), observed effects 
from geometrical nonidealities in this study have important implica- 
tions. It shows that suspended 2D structures appearing to be 'imper- 
fect' should not be deemed unworthy, instead they make very 
interesting devices that can offer new functions (e.g., enabled by 
the free edges). Collectively, these 2D structures with 'nonideal' 



SCIENTIFIC REPORTS | 4:3919 | DOI: 1 0.1 038/srep0391 9 



5 



appearances or broken symmetries form new families of resonators 
with versatile dispersion relations (of multimode frequency versus 
geometry). It also demonstrates that geometry can be powerful for 
tuning the device performance, such as resonance frequencies, mode 
shapes, mode splitting, and frequency spacing. It is therefore intri- 
guing to engineer devices with 'customized' resonance behavior, such 
as specified frequency ratios between the different resonance modes, 
by tailoring device geometric parameters. This can be a useful 
tool, for example, in designing multimode and nonlinear resonant 
components in micromechanical circuits, where the ability of cre- 
ating arbitrary ratios between resonance frequencies of the different 
modes is desirable 54 . The capability of controlling mode shapes can 
also be useful in designing resonators with position- specific surface 
adsorption sensing and manipulation capabilities 49,51 . We conclude 
that visual appearance of MoS 2 and other 2D devices should not be a 
criterion forjudging the devices; embracing and exploiting structural 
nonidealities and broken symmetries can benefit future device inno- 
vations. Looking forward, devices with these interesting structural 
nonidealities will not be limited to random yield in hand- craft exfo- 
liation on small chips, but can be intentionally designed - such 
opportunities may soon be enabled by advancements in high-pre- 
cision deterministic transfer methods 7 " 11 and large-area synthesis of 
MoS 2 with improving control (e.g., of size, shape, crystalline 
orientation) 55 " 58 . 

Methods 

Device fabrication. The devices are fabricated using the mechanical exfoliation 
method. First, we perform optical lithography followed by buffered oxide etch (BOE) 
to pattern arrays of microscale trenches of various sizes and shapes on a 290 nm-thick 
silicon dioxide (Si0 2 ) layer thermally grown on a silicon (Si) wafer. The BOE etch 
time is controlled such that the entire depth (290 nm) of Si0 2 is removed, exposing 
the flat Si surface underneath. We then exfoliate MoS 2 flakes onto the structured 
wafer. Flakes covering the microtrenches form drumhead -structured 
nanomechanical resonators. 

Mechanical exfoliation of MoS 2 typically results in small flakes (i.e., often smaller 
compared to exfoliated graphene flakes on the same substrate using exactly the same 
procedure). When the flakes are small (with respect to the microtrenches), they are 
likely to form devices where the microtrenches are not fully covered by the MoS 2 
nanosheets. Incomplete drumheads with partially-clamped boundaries occur, cre- 
ating incompletely- clamped drumhead nanomechanical resonators. Supplementary 
Fig. SI shows a schematic of the device fabrication process and the resulting device 
geometries. 

Interferometric motion detection system. The mechanical resonances of the 
nonideally or imperfectly structured MoS 2 resonators are measured using a custom- 
built system with ultrasensitive optical readout techniques 21 ' 22 . The schematic of the 
measurement system is shown in Supplementary Fig. S2. 

In our custom-built interferometry system, the laser beam from a 633 nm He-Ne 
laser first goes through a beam expander. It is then focused by a 50 X microscope 
objective lens with optical window correction (NA = 0.5) onto the surface of the 
MoS 2 nanosheet (the spot size is — 1 urn). The interference between the light waves 
reflected from the suspended MoS 2 nanosheet and the underneath Si surface deter- 
mines the intensity of the total reflected light. The undriven thermomechanical 
motions (i.e., Brownian motions, as dictated by the fluctuation-dissipation theorem) 
of the MoS 2 nanosheet modulate the interference signal. The interferometric signal is 
detected by a low-noise photodetector (PD), and then sent to a spectrum analyzer. 
The high sensitivity of the measurement system allows us to detect the minimal 
vibrations of the MoS 2 nanosheet arising from its fundamental thermomechanical 
(Brownian) motions and thus find its resonances without any external driving. We 
also adjust our on-device laser power to be <700 uW such that the heating effect is 
negligible (no heating- induced frequency change within the power range used), while 
maintaining a good signal-to-noise ratio 22 . All our measurements are conducted in 
moderate vacuum (~6 mTorr) at room temperature. 
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